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The theoretical explanations about the “pa puzzle” in charmonium decays are reviewed exten- 
sively, and the comparison of theoretical predications with experimental data is made whenever 
possible. Three methods to estimate the ratio of the branching fractions of J/y and w’ decays are 
also discussed. It is pointed out that in order to understand the pa puzzle, and the dynamics of 
charmonium decays, systematic studies should be made in theory, phenomenology and experiment 


aspects. 
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I. INTRODUCTION 


Crisply defined experimental puzzles in high-energy 
physics always have the prospect of leading to new discov- 
eries. One prominent example is the 6-7 puzzle of 1956 
which led to the parity revolution. Therefore puzzles in 
physics often draw considerable attention to theorists. 
The ratios of hadronic decays of the w(3686) (shortened 
as w’) to the same decays of the J/w is such a puzzle 
which has been studied substantially since 1983. 

Since the OZI suppressed decays of J/y and 7’ to 
hadrons are via three gluons or a photon, in either case, 
the perturbative QCD (pQCD) provides a relation [1] 


Byrn — By sete- 


Qn = x 12.7% . (1) 


Bjjy—n 7 Byjpete- 


This relation is referred to as the “12% rule” which is 
expected to be held to a reasonable good degree for both 
inclusive and exclusive decays. The so-called “pa puzzle” 
is that the prediction by Eq. (1) is severely violated in 
the pa and several other decay channels. The first evi- 
dence for this effect was found by Mark-II Collaboration 
in 1983 [2]. From then on many theoretical explanations 
have been put forth to decipher this puzzle. 

With the recent experiment results from BESII 
and CLEOc about J/wW and w’ two-body decays, 
such as vector-pseudoscalar (VP), vector-tensor (VT), 
pseudoscalar-pseudoscalar (PP), and baryon-antibaryon 
(BB) modes, and about multi-body decays at the J/1, 
the 7’ or even at the ~(3770) (shortened as %”) [3]-[23], a 
variety of solutions proposed for the puzzle can be tested 
at the level of higher accuracy. In this treatise, we survey 
the theoretical works on the pr puzzle and compare them 
with the available experimental data. From the theoret- 
ical point of view, since the Q-value is smaller than 12% 
for pr, it may be caused either by enhanced or suppressed 
J/w decay rate. Another possibility is by both. So we 
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classify the relevant theoretical speculations into three 
categories: 


1. J/w-enhancement hypothesis, which attributes the 
small Q-value to the enhanced branching fraction 
of J/w decays. 


2. ~’-suppress hypothesis, which attributes the small 
Q-value to the suppressed branching fraction of Y’ 
decays. 


3. Other hypotheses, which are not included in the 
above two categories. 


In the following content, first reviewed are the theoreti- 
cal works on pr puzzle, and the predictions from them are 
compared with the newly available experimental results; 
then expounded are three methods to estimate the ratio 
between w’ and J/w decays into the same final states; 
after that some comments are made on the implications 
from the review on pr puzzle; at last there is a short 
summary. 


II. REVIEW OF THEORETICAL WORKS ON pr 
PUZZLE 


A. J/i-enhancement Theory 


In the earlier days of the pa puzzle, it was noticed 
that the decay of 1~~ charmonium into pr final state vio- 
lates the Hadronic Helicity Conservation (HHC) theorem 
(see below for expound) [24], and so such decay should 
be suppressed. Therefore people think there must be 
some mechanism which leads to the great enhancement 
for J/w — pr decays. The two schemes presented in this 
section were proposed following this line of reasoning. 


1. J/tb-glueball Admixture Scheme 


The idea of J/y decays via a glueball was proposed 
by Freund and Nambu [25] (FN hereafter) soon after 
the discovery of J/~ particle to explain the breaking 


of Okubo-Zweig-lizuka (OZI) rule [26]. In such mech- 
anism, the breaking results from the mixing of the w, 
@, and J/w mesons with an SU(4)-singlet vector meson 
O. They found that such an O meson should lie in the 
1.4-1.8 GeV/c? mass range with the width greater than 
40 MeV /c?, and it should decay copiously into pr, K*K 
while exhibiting severe suppression of decays into KK, 
ete” and utu modes. These authors presented several 
quantitative predications for experimental search. Two 
of them are 


T 2 i 
R= — 0.0115 — 0.087 , 
Lgpr 
T —_ 
Ry = ee <8 x 10%. 
J/b—pr 


With the current available data and using the three pi- 
ons final state as a substitute for pr in both ¢ [27] and 
J/w [3, 4] decays, we obtain the first ratio R ~ 0.003, 
which is almost one order of magnitude smaller than 
the predication. For the second ratio, by virtue of 
PDG [27] value for K+K- and new experiment result for 
KK? [6], it is estimated that B(J/y > KK) ~ 1074, 
together with the results for pr [3, 4], we have Ro ~ 107? 
which is much larger than the predication. 

The first attempt to explain the pr puzzle in terms of a 
glueball near J/¢ was proposed by Hou and Soni [28] (HS 
hereafter). They attributed the enhancement of J/y —> 
K*K and J/ — pr decay modes to the mixing of the Y 
with a JPC = 177 vector gluonium, also designated by 
O. The differences between FN’s and HS’s pictures lie in 
the following aspects: 


e Based on the potential model applied to glueball, 
the mass of a low-lying three-glue state is estimated 
to be around 2.4 GeV/c? [29], rather than 1.4 to 1.8 
GeV/c? in Ref. [25]. 


e The mixing of O with w’ is taken into account, 
which has been ignored in previous work. 


Since the gauge coupling constant in QCD is mo- 
mentum dependent, the mixing parameter is taken 
to be a function of the invariant mass q?, which de- 
creases rather sharply with the increase in q?. Such 
propagator effect gives rise to a more suppression 
on the decay rates of Y’ relative to J/w for decays 
to pr and K*K channels. 


By virtue of their assumption, HS suggested a search 
for the vector gluonium state in certain hadronic decays 
of the Y’, such as Y’ > nr +X, 7+X, n +X, where X 
decays into VP final states [28]. 

Based on HS’s idea, Brodsky, Lepage, and Tuan [30] 
(BLT) refined the glueball hypothesis for the pr puzzle. 
They assumed the general validity of the pQCD theo- 
rem that the total hadron helicity is conserved in high- 
momentum-transfer exclusive processes, in which case 


the decays to pr and K*K are forbidden for both the J/w 
and w’. This pQCD theorem is often referred to as the 
rule of Hadronic Helicity Conservation (HHC) [24], which 
is based on the assumption of the short-range “pointlike” 
interactions among the constituent quarks throughout. 
For instance, J/y(ct) — 3g has a short range ~ 1/me 
associated with the short time scale of interaction. Nev- 
ertheless, if subsequently the three gluons were to res- 
onate forming an intermediate gluonium state O which 
has a large transverse size covering an extended time pe- 
riod, then HHC would become invalid. In essence the HS 
model takes over in this latter stage. 

Final states h which proceed only through the inter- 
mediate gluonium state satisfy the ratio 


G= Bib! — ete-) (Myjy — Mo)? +T3/4 2 
”— B(JJp > ete-) (My — Mo)? +T3/4 ` 


NN 


The Qn is small if the mass of O is close to the mass 
of J/w. The experimental limits at that time [2, 30, 31] 
imply that the O mass is within 80 MeV/c? of the mass of 
J/w and its total width is less than 160 MeV/c?. Brodsky 
et al. recommended a direct way to search for O, that 
is to scan the ete~ — VP cross section across the J/w 
resonance. 

Another related work by Chan and Hou [32] studied 
the mixing angle 004 and the mixing amplitude foy of 
the J/w and vector glueball O based on the framework 
of potential models of heavy quarks and constituent glu- 
ons. They obtain |tan@oy| = 0.015 and foy(ms,) = 


0.008 GeV’. 

On the experimental part, BES has searched for this 
hypothetical particle in a pr scan across the J/w mass 
region in e*e~ annihilations as well as in the decays of 
yw’ — ntO, O — pr, and found no evidence for its ex- 
istence [8, 33]. The data constrains the mass and width 
of the O to the range | Mo — Myyy, |< 80 MeV/c? and 
4 <To < 50 MeV/c? [34]. Although the absence of dis- 
tortion in BES energy scan of J/y~ — pr does not rule 
out Mo ~ M74, it puts a lower bound to ro. However, 
as indicated in Ref. [36], the experimentally constrained 
mass is several hundred MeV/c? lower than the mass of 
the lightest vector glueball calculated in lattice simula- 
tions of QCD without dynamical quarks [37]. 

Recently, the experimental data from BES and CLEOc 
turned out to be unfavorable to this glueball hypothesis. 
Among them is the observed large branching fractions 
of the isospin-violating VP mode 7’ — wr? [9, 10, 16]. 
This contradicts with the assertion that the pattern of 
suppression depends on the spin-parity of the final state 
mesons. In addition, according to BLT’s analysis, one 
obtains the relation [38] 


B(J/ġ > wr) 


— ~ < 0.0037 
BOJ > Pr) ~ 


which is much smaller than the PDG06 value 0.08 [27]. 
Another experimental result which is unfavorable to such 


hypothesis is the suppression of 7’ decays into vector- 
tensor (VT) final states [18, 19]. Since hadronic VT de- 
cays, unlike the VP decays, conserve HHC, some other 
mechanisms must be responsible for this suppression in 
the model. Furthermore, it has been argued that the O 
may also explain the decay of J/ẹ into fo (named pre- 
viously S*) but not to pag(980) (named previously ô), 
since the O mixes with the ¢@ and enhances a mode that 
would be otherwise suppressed [30]. However, the decay 
w' — fo [34] is not suppressed as found by experiments, 
which implies the absence of anomalous enhancement 
in J/y) — fo, thus contradicts with this explanation. 
Anselmino et al. extended the idea of J/W-O mixing to 
the case of ne + VV and pp [39]. They suggested that 
the enhancement of these decays can be attributed to the 
presence of a trigluonium pseudoscalar state with a mass 
close to the 7. mass. So far there is no experimental 
evidence for such a state. 

In fact, as pointed out in Ref. [40], this glueball ex- 
planation has some unanswered questions: (i) Why only 
the pr and K*K channels are affected but not O — 57 
etc. viz., one must assume in an ad hoc way that the O 
couples predominantly to pr and K*K; and (ii) if such 
a narrow heavy 1~~ gluonium state exists, why have the 
narrow 0+*, 2++ states not been found, which should be 
lighter and easier to search? 


2. Intrinsic-charm-component Scheme 


Brodsky and Karliner (BK) put forth an explana- 
tion for the puzzle based on the existence of the in- 
trinsic charm |ggc¢) Fock components of the light vector 
mesons [41]. They noticed the fact that quantum fluctua- 
tions in a QCD bound state wave function will inevitably 
produce Fock states containing heavy quark pairs. The 
intrinsic heavy quark pairs are multiconnected to the va- 
lence quarks of the light hadrons, and the wave functions 
describing these configurations will have maximal ampli- 
tude at minimal off-shellness and minimal invariant mass. 
In the case of the p meson, with the consideration of the 
light-cone Fock representation: 


p™ = Yf-lud) + Yf- ludce) pats 


Here we expect the wave function of the cc quarks to be 
in an S-wave configuration with no nodes in its radial de- 
pendence, in order to minimize the kinetic energy of the 
charm quarks and thus also minimize the total invariant 
mass. 

The presence of the |udcé) Fock state in the p allows 
the J/y — pr decay to proceed through rearrangement 
of the incoming and outgoing quark lines; in fact, the 
|udcc) Fock state wave function has a good overlap with 
the radial and spin |c@) and |ud) wave functions of the 
J/w and pion. On the other hand, the overlap with the 
w’ will be suppressed, since the radial wave function of 
the n = 2 quarkonium state is orthogonal to the nodeless 
cT in the |udc@) state of the pr. Similarly, the usce) Fock 


component of the K* favors the J/wK configuration, al- 
lowing the J/y — K*K decay to also proceed by quark 
line rearrangement, rather than c¢ annihilation. 

These authors also suggested comparing branching 
fractions for the 7. and 7, decays as clues to the im- 
portance of 7. intrinsic charm excitations in the wave- 
functions of light hadrons. 


B. y’-suppress Theory 


The hypothesis of the existence of a glueball to explain 
the pr has been questioned soon after it was proposed. 
In addition, it is also pointed out [42] that the helicity 
suppression is not a strong constraint in the charmonium 
energy scale. Under such case, one comes naturally to 
the idea that it is not J/W — pr which is enhanced, but 
rather ~’ — pr which is suppressed. Seven explanations 
or models collected in this section are put forth along this 
line. 


1. Sequential-fragmentation Model 


Karl and Roberts have suggested explain the pr puz- 
zle based on the mechanism of sequential quark pair cre- 
ation [43]. The idea is that the quark-antiquark pairs 
are produced sequentially, as a result the amplitude to 
produce two mesons in their ground state is an oscilla- 
tory function of the total energy of the system. They ar- 
gue that the oscillatory fragmentation probability could 
have a minimum near the mass of w’, which provides an 
explanation for the suppressed wy’ decay. Even though 
their evaluations could generally accommodate the data 
for decays of J/y and y’ to pr and K*K, it runs into 
difficulties when it is extrapolated to T decays. Accord- 
ing to their calculation, the oscillations of probability 
amplitude are damped out in the region of the T res- 
onances, so the pr channel is present in the decay of all 
Y, Y’, T”, --- resonances with a common rate. This 
leads to a prediction [(Y — pr) = 0.05 keV, or equiv- 
alently B(T — pr) = 9.4 x 1074, which is above the 
current upper limit B(Y — pr) < 2 x 1074 [27]. More- 
over, their calculation seems hard to explain the large 
branching fraction for ¢ decays to pr [27] due to the fact 
that their fragmentation probability tends to zero as the 
mass of the pr decaying system approaches 1GeV/c?. 

In a further analysis [44], Karl and Tuan pointed out 
that if a suppression is observed in three-meson channels 
the explanation based on sequential pair creation would 
be undermined. Recently such a suppressed channel, viz. 
@KK, is found by CLEOc [21]. 


2. Exponential-form-factor Model 


Guided by suppressed ratios of y’ to J/wW decays to 
two-body hadronic modes, Chaichian and Törnqvist sug- 


gested [40] that the hadronic form factors fall exponen- 
tially as described by the overlap of wave functions within 
a nonrelativistic quark model. This behavior explains 
the drastically suppressed two-body decay rates of the 
w’ compared with those of the J/w. Recently, the re- 
port on observation of a number of VP channels in y’ 
decays [9, 10, 16] such as wy’, dn’, pn’ has proved that 
the predicted decay fractions are overestimated. More- 
over, the branching fraction for wr? [27], is well below 
the prediction by this model which is 1.04 x 1074. 

Another problem of the model is that it does not single 
out just the VP channel, the other channels, for example 
VT channel, are also estimated to have small branching 
fractions which are not compatible with the BES mea- 
sured results [19]. 


3. Generalized Hindered M1 Transition Model 


A so-called generalized hindered M1 transition model 
is proposed by Pinsky as a solution for the puzzle [45]. 
It is argued that because J/y — yn is an allowed M1 
transition while y’ — yn is hindered (in the nonrela- 
tivistic limit), using the vector-dominance model to relate 
yy! — y to Y’ — wr! one could find the coupling Gyn, 
is much smaller than Gy, ,),, and then by analogy, the 
coupling Gypr would be much smaller than Gwpr. Here 
Guwpr can be extracted from data by virtue of the anal- 
ysis using the vector-dominance model and a standard 
parameterization of OZI process [46]. Then assuming 
w’ — pr to proceed via wW’-w’ mixing, while J/Y — pr 
via ~-w mixing, one would find that Yy’ — pr is much 
more severely suppressed than J/y — pr. The similar 
estimation could be preformed for K*K and other VP 
final state, and one can expect a suppressed Q: 


B(Y > VP) B Troal) (every Fy 
B(y > VP) © Pial) \Gvve / Fy 


= 0.06% , 


(3) 
where Fy /Fy = 0.3, Gurpr/Gupr = 0.066 according 
to Ref. [45]. This Q is much smaller than the present 
experimental results [9, 10, 16]. 

Moreover, in this model, the coupling Gw'wfa for w” > 
wf should not be suppressed because by analogy the 
coupling Gyyy.. is not small due to the fact that the E1 
transition w’ —> YXe2 is not hindered [47]. Therefore via 
y’-w’ mixing the Y — w’ — w fz decay is expected to be 
not suppressed, which contradicts the BES result [19]. 


4.  Higher-Fock-state Scheme 


Chen and Braaten (CB) proposed an explanation [36] 
for the pr puzzle, arguing that the decay J/w — pr is 
dominated by a Fock state in which the cc pair is in a 
color-octet 3S4 state which decays via c¢ — qq, while the 
suppression of this decay mode for the 7’ is attributed to 
a dynamical effect due to the small energy gap between 


4 


the mass of the y’ and the DD threshold. Using the 
BES data on the branching fractions into pr and K*K 
as input, they predicted the branching fractions for many 
other VP decay modes of the w’, as listed in Table I, from 
which we see most measured values falling in the scopes 
of predictions, but we also note for wa mode, the devi- 
ation from the prediction is obvious. Here it should be 
noticed that the values adduced in Table I are calculated 
on the strength of the measured branching fractions from 
earlier experiment, the new measurements on the branch- 
ing fractions for pr and K°R® + c.c. from BES [11, 14] 
and CLEOc [16] may have impact on the corresponding 
evaluations. 


TABLE I: Predictions and measurements for Qvp in unit of 
1% for all VP final states. The value for pr and ROR tae: 
from Ref. [48] were used as input. The theoretical parameter 
x = 0.64 is due to results of Ref. [49] and the experimental 
results come from Ref. [9-11, 50]. 


VP x = 0.64 Exp. 
pr 0—0.25 0.13 Æ 0.03 
K*° K? + c.c. 1.2 — 3.0 3.2 0.08 
K*+K~ +c.c. 0— 0.36 0.59+9:27 


wn 0-16 <20 
wn 12—55 19433 
on 0.4-3.0 5.11.9 
on 0.5-2.2 9.44.8 
pn 14—22 9.2435 
pry 12—20 17.841? 9 
wT 11—17 4.4tT 8 


Besides the predictions in Table I, CB’s proposal also 
has implications for the angular distributions for two- 
body decay modes. In general, the angular distribution 
has the form 1 + acos?6, with —1 < a < +1. CB’s 
conclusion implies that the parameter a for any two-body 
decay of the 7’ should be less than or equal to that for 
the corresponding J/w decays. But this needs further 
supports from analyses based on large data sample in 
the future. 


5. Survival-chamonia-amplitude Explanation 


A model put forward by Gérard and Weyers enter- 
tains the assumption that the three-gluon annihilation 
amplitude and the QED amplitude add incoherently in 
all channels in J/w decays into light hadrons, while in the 
case of yy’ decays the dominant QCD annihilation ampli- 
tude is not into three gluons, but into a specific config- 
uration of five gluons [51]. More precisely, they suggest 
that the strong annihilation of the ~’ into light hadrons is 
a two-step process: in the first step the Y’ goes into two 
gluons in a 0** or 07+ state and an off-shell h.(3526); 
in the second step the off-shell he annihilates into three 
gluons to produce light hadrons. Their argument implies: 


(a) to leading order there is no strong decay amplitude 
for the processes y’ — pr and Yy’ > K*K; (b) the 12 % 
rule should hold for hadronic processes which take place 
via the QED amplitude only. As far as the second im- 
plication is concerned, the present data give different ra- 
tios between Y’ and J/y decay for wn? and mtr final 
states, both of which are electromagnetic processes. Here 
the event form factor effect is taken into account [52], the 
difference between two kinds of processes is still obvious. 
Besides the explanation for pr puzzle, this model predicts 
a sizable Y’ — (ntn or 7) hy(1170) branching fraction. 

In a recent paper [53], Artoisenet, Gérard and Wey- 
ers (AGW) update and sharpen the above idea which 
leads to a somewhat unconventional point of view: all 
non-electromagnetic hadronic decays of the y’ goes via 
a transition amplitude which contain a cé pair. AGW 
provide two patterns for these two-step decays, the first 
is 


wy’ — (2NPg) + (39) . (4) 


The physics picture is as follows: the excited cc pair in 
the 7’ does not annihilate directly. Instead, it spits out 
two non-perturbative gluons (2NPg) and survives in a 
lower cœ configuration (1~~ or 1~+) which then eventu- 
ally annihilate into 3g. The decays y’ — (27)J/y and 
w’ + nJ/w follow this pattern. The second pattern is 


wy’ — (3NPg) + (2g) , (5) 


where the lower c¢ configuration (07* or 0**) annihilates 
into 2g. The only on-shell channel for this type of decays 
is Y’ — (37), whose branching fraction is estimated as 
(1 — 2)% level. Anyway, the recent measurement from 
CLEOc [20] provides the upper limit which is one order 
of magnitude below this theoretical prediction. Further- 
more, the substitution of one photon for one gluon in 
Eqs. (4) and (5) allows 


wy’ — (2NPg) + (29) +7. (6) 


This pattern corresponds to on-shell radiative decays 
such as Y —> (ntr )ney and Y —> ney, which could 
be larger than the observed Y’ — ney mode. 

Besides the above predications, AGW also estimate 


B(w' > bin) = (1.3 0.3) x 1073 , (7) 
Bi! > hın?) = (1.9 0.4) x 107% , (8) 
BUJ/b > bin) ~ B(Y > bin) x1% . (9) 


All these are to be tested by experiments. 


6. Nonvalence Component Explanation 


Since the w’ is a highly excited state and close to the 
DD threshold, it is suggested [54] that unlike the J/1, 
the w’ may be an admixture of large nonvalence com- 
ponents in its wave function. The so-called nonvalence 


component indicates the additional gluon or the light 
quark-antiquark pair (or as in Ref. [54], the cég compo- 
nent and a quasi-molecular DD state), which makes 7’ 
decays prominently distinctive from those of J/w. The 
nonvalence component of the J/w is expected to be less 
significant because it is the lowest state. In a following 
paper [42], Chernyak uses this picture to explain the pr 
puzzle. He suggested that the valence and nonvalence 
strong contributions interfere destructively in pa chan- 
nel and consequently cancel to a large extent in the total 
p! — pr strong amplitude, while the role of nonvalence 
contributions is much less significant in J/w — pr. From 
this viewpoint, there is no deep reason for the experimen- 
tally observed very strong suppression of Yy’ — pr, this 
is the result of a casual cancellation. 

Chernyak also tries to use the above idea to explain 
qualitatively other decay modes, such as VT, AP, PP, 
VV and BB decay. However, such an idea remains a 
pure speculation, and no concrete calculations have been 
performed up to now. 


7 S-D Wave Mixing Scheme 


The 4%” is viewed as a D-wave charmomium state. Al- 
though it is primarily cc(1°Dj), its leptonic width in- 
dicates a contribution from mixing with S-wave states, 
mainly the nearby 7(23,$1). This leads Rosner to propose 
that the small pr branching fraction in w’ decay is due 
to the cancellation of the 2S and 1D matrix elements. 
By virtue of his scheme 


(om |’) = (pr|23S1) cos@ — (pr|13Dj) sin @ , 
lomld’) = (pm|23S,) sin 8 + (pm|13D;) cos9, O 
where @ is the mixing angle between pure 7(2°.$,) and 
w(13D,) states [55] and is fitted from the leptonic widths 
of the ~” and the y’ to be (12+ 2)° [56], which is consis- 
tent with the coupled channel estimates [57, 58] and with 
the ratio of w’ and y” partial widths to J/aata— [59]. 
If the mixing and coupling of the y’ and y” lead to com- 
plete cancellation of Yy’ — pr decay ((p7|~") = 0), the 
missing pr decay mode of the y’ shows up instead as de- 
cay mode of the y”, enhanced by the factor 1/ sin? 0, the 
concrete estimation shows that [56] 


Byr—pr = (4.1 + 1.4) x 1074 . (11) 


To test this scenario by the data collected at the w” in 
ete” experiments, it has been pointed out [60, 61] that 
the continuum contribution must be considered carefully. 
Specifically speaking, by Rosner’s estimation, the Born 
order cross section for Y” — pr is 


ogo n = (4.8 + 1.9) pb, (12) 


— pt 


which is comparable in magnitude to that of the contin- 
uum cross section, viz. 


of. =4.4 pb. (13) 


ete-—pt 


So what is observed is the coherent sum of the two am- 
plitudes. Based on the analysis of experimental data, 
it has been suggested that there be a universal phase 
between strong and electromagnetic amplitudes in char- 
monium decays. With this assumption, the strong decay 
amplitude interferes with the continuum amplitude ei- 
ther maximum destructively, e.g. for pr, wy and K*+ K- 
or maximum constructively, e.g. for K*°K°. The de- 
structive interference leads to the phenomena that the 
measured cross section on top of the resonance could be 
smaller than that off the resonance. Recent experimental 
results [13, 22] on pr, wn and K** K- have demonstrated 
such interference pattern. This provides support to the 
prediction of Eq.(11). However, to uniquely determine 
By on, current available experimental data are insuf- 
ficient, the Y” resonance must be scanned [62]. So the 
quantitative test of Rosner’s scenario remains to be the 
task of the future experiments on t-Charm factories. 

In the consequent study [63], Wang, Yuan and Mo 
(WYM) extend the S-D wave-mixing scenario to PP fi- 
nal state, and give a constraint for Y” > K&K? decay, 


0.12 + 0.07 < 10° x Bi” > KĮK?) < 3.8411, (14) 


which is compatible with both the BESH result B(w” — 
KLK?) <2.1x 1074 at 90% C.L. [12] and the CLEOc re- 
sult B(Y” > KLK?) < 117x107 at 90% C.L. [23]. Ex- 
trapolating fore-study to all charmless decays [64], WYM 
found that this scenario could accommodate large non- 
DD decay of the y”. Although the recent experimental 
measurement from CLEOc favors the zero non-DD cross 
section at Ww” [65], the larger error could not rule out the 
existence of non-DD branching fraction at 10% level. 


C. Other Explanations 


Besides the models in two previous sections, more spec- 
ulations involving the pa puzzle are collected in this sec- 
tion. 


1. Final State Interaction Scheme 


Li, Bugg and Zou [66] (LBZ) pointed out that the fi- 
nal state interactions (FSI) in J/w and y’ decays give 
rise to effects which are of the same order as the tree 
level amplitudes, they argued that J/w — pr is strongly 
enhanced by the agp loop diagram, while the direct tree- 
production for pz may be suppressed by the HHC mech- 
anism [24]. The contribution of the ap loop diagram is 
much smaller than that of agp loop for the J/w — pr, but 
they have similar strength for the Yy’ — pr and may can- 
cel each other by interfering destructively. The similar 
apparent suppression for y’ — K*K and fow may also 
be explained by the K “Kaa and bır loop, respectively. 
Therefore, LBZ expected that FSI may provide a coher- 
ent explanation for all the observed suppressed modes of 


wy’ decays. However, as remarked in Ref. [67], this inter- 
ference model does appear to have more assumptions that 
predictions. The only qualitative prediction by LBZ is to 
check whether the ajp and Kj K production rates are 
large for the y’. So far, no such measurements have been 
reported. Nevertheless, useful information on ap and 
iy K could be obtained from the BES published data 
as shown in Refs. [17] and [19]. The lack of evidence 
within the invariant mass distribution plots (see Fig.3 
and Fig.5 of Ref. [17] or Ref. [19]) that the pr recoiled 
against a p for events of Y’ — p°p*nt and that + KF 
recoiled against a K*° for events of Y > ntr- KK- 
suggests that they are unlikely to be the favored modes 
in w’ decays. 


2. Large Phase Scheme 


Suzuki gave another explanation based on FSI for J/w 
decays [38]. He performed a detailed amplitude analy- 
sis for J/W — 1707 decay to test whether or not the 
short-distance FSI dominates over the long-distance FSI 
in the J/w decay. His result indicates that there is a 
large phase between three-gluon and one-photon ampli- 
tudes. Since the large phase cannot be produced with the 
perturbative QCD interaction, the source of it must be 
in the long-distance part of strong interaction, namely, 
rescattering among hadrons in their inelastic energy re- 
gion. Suzuki then performed the similar analysis for 
J/w — 07-07 decay, and obtained the same conclusion 
of the large phase [68]. His analysis also shows that the 
exclusive decay rate at the J/w is in line with that of 
the inclusive decay. This fact leads him to believe that 
the origin of the relative suppression of 7’ — 1707 to 
J/w — 1707 is not in the J/w but in the y”. 

As to this large phase, in Ref. [51], Gérard argued that 
this phase follows from the orthogonality of three-gluon 
and one-photon virtual processes. As a matter of fact, the 
conclusion of a large phase close to 90° between three- 
gluon and one-photon processes holds true for all two- 
body decays of J/~, such as 1707 (90°) [69], 1707 (1064 
10)° [70, 71], 1717 (138 +37)° [68, 71, 72], 0707 (89.64 
9.9)° [68, 71, 72] and NN (89+ 15)° [71, 73]. 

Nevertheless, as the Y’ decays were analyzed, the ex- 
perimental data at first seemed to favor a phase close to 
180° [69], in contrary to the expectation that the decay 
of J/w and w’ should not be much different. However, 
as pointed out by Wang et al. that the previous pub- 
lished data did not take the continuum one-photon pro- 
cess into account. Their reanalysis of the experimental 
data shows that the phase with value around —90° could 
fit y’ > 1707 data [61] and +90° could fit Yy’ — 0707 
data [74]. The latter is confirmed by more recent re- 
sults by CLEO [5]. Furthermore, this large phase also 
shows in the OZI suppressed decay modes of the y”. 
In many decays modes of the 7”, the strong decay am- 
plitudes have comparable strength as the non-resonance 
continuum amplitude, the large phase around —90° leads 


to destructive or constructive interference. Due to the 
destructive interference, the observed cross sections of 
some modes at the peak of the 7” are smaller than the 
cross section measured off-resonance [60]. This is demon- 
strated by the data from CLEOc [22]. 

If the large phase between the three-gluon and one- 
photon virtual processes is universal not only in J/ de- 
cays, but also in y’ as well as all charmonium or perhaps 
all quarkonium decays, then how this phase is related to 
the difference between the decays of J/y and Y’ remains 
to be an unanswered question. 


3. Mass Reduction Explanation 


In the study [75] of radiative decays of 1~~ quarko- 
nium into 7 and 7’, Ma presented a QCD-factorization 
approach, with which he obtained theoretical predictions 
in consistency with CLEOc measurement. The largest 
possible uncertainties in analysis are from the relativistic 
corrections for the value of the charm quark mass. Ma ar- 
gued that the effect of these uncertainties can be reduced 
by using quarkonium masses instead of using quark mass. 
As an example of such reduction, he provided a modified 
relation to the original 12% rule 


Q- BU = er) _ Mi )y BJ / > ete-) 
"BC! > pr) M, BY! > ete) 
= (3.6 + 0.6)% . 


However, this value is much larger than the experimental 
result given in Table I. 


4. Vector-meson-mixing Model 


Intending to give a comprehensive description of 
J/w two-body decays, Clavelli and Intemann (CI) pro- 
posed a vector-meson-mixing model in which the vector 
mesons (p, w ,¢@ ,J/w) are regarded as being admix- 
ture of light-quark-antiquark state and charmed-quark- 
antiquark state [76]. The coupling of the J/7 to any 
state of light quarks is then related to the corresponding 
coupling of the p, w, and ¢ to the same state. With few 
experiment inputs to determine the mixing parameters, 
CI calculated VP, PP, and BB decay rates for the J/y as 
a function of the pseudoscalar mixing angle. Most of the 
predictions agree with the experiment results at the order 
of magnitude level, but discrepancy is obvious for some 
channels, such as KK? final state [6]. CI also extended 
their model to the hadronic decays of the y’. Neverthe- 
less, their evaluations for B(J/y —> wn?) = 3 x 1075 and 
B(W' — wr?) = 3 x 107% contradict with the present 
results (4.5 + 0.5) x 1074 and (2.1 + 0.6) x 1074 [27], 
respectively. 

Starting from effective Lagrangian whereby nonet- 
symmetry breaking and pseudoscalar-meson mixing can 


be studied, Haber and Perrier parametrized the decay 
modes of J/y — PP (also for J/y => VV or qe > 
VP), J/y — VP (also for J/W — VT or ne > VV), 
J/w — PPP (also for J/y = VVP or ne > PPV), and 
Ne > PPP (also for J/y — PPV or ne > VV P) [77]. 
The experimental data are used to determine these phe- 
nomenological parameters. In a later work, Seiden, 
Sadrozinski and Haber took the doubly Okubo-Zweig- 
lizuka suppression (DOZI) effect into consideration, and 
presented a more general parameterization of amplitudes 
for J/y) — PP decays [78]. With this form, one could 
easily derive the relative decay strength between different 
final states. However, it has also been noticed that un- 
der the most general circumstances, symmetry breaking 
introduces too many parameters for a general analysis to 
be useful. In order to reduce the number of new parame- 
ters considerably and make the analysis manageable, only 
special cases could be considered. 


TABLE II: Comparison of predictions [79] and measure- 
ments [27] for the branching ratios (10~°) for y’ decays. The 
upper limits are presented at 90% C.L. 


VP Prediction Measurement 
pr 1.3 3.2+1.2 
KR + c.c. 5.1 10.9 + 2.0 
K**K7>+ece 13 i aba 
wn 1:2 < 1.1 
wn’ 6.3 20°05 
on 1.6 25h e 
on! 4.6 3.11.6 
pn 21 2.2+0.6 
pr! 1:3 157s 
wr 3.8 2.1+0.6 
or 0.01 < 0.40 


A similar parameterization with mixing feature of the 
strong interaction mechanism was proposed by Feld- 
mann and Kroll (FK) [79] for the hadron-helicity non- 
conserving J/~ and y’ decays, but with a different inter- 
pretation from those put forth in Refs. [36, 49, 67, 78]. 
FK assumed that with a small probability, the charmo- 
nium possesses Fock components built from light quarks 
only. Through these Fock components the charmonium 
state decays by a soft mechanism which is modeled by 
J/w-w-d mixing and subsequent w (or ¢) decays into the 
VP state. In absence of the leading-twist perturbative 
QCD contribution, the dominant mechanism is supple- 
mented by the electromagnetic decay contribution and 
DOZI violating contribution. FK argued that this mech- 
anism can probe the charmonium wave function at all 
quark-antiquark separations and feels the difference be- 
tween a 1S and a 2S radial wave functions. The node 
in the latter is supposed to lead to a strong suppression 
of the mechanism in ~’ decays. With a few parameters 
adjusted to the experimental data, FK obtained a nu- 
merical description of the branching fractions for many 
VP decay modes of the J/y and Y’, which agree with the 


measured branching fractions at the order of magnitude 
level, as shown in Table II. Moreover, FK has extended 
their mixing approach to the ne — VV decays and ob- 
tained a reasonable description of the branching fractions 
for these decays while the 7/, + VV decays are expected 
to be strongly suppressed. 


D. Remarks 


From the brief retrospect of pa puzzle history, we 
could notice that in the early stage, theorists concen- 
trate on peculiar mechanism for special channels like pa 
and K*K, such as Hou and Soni’s explanation. With 
the development of the theory and the progress of the 
experiment, theorists attempt to provide a more general 
scheme for charmonium decays, such as the work by Feld- 
mann and Kroll. In fact, the charmonium decay is an 
interconnected system as a whole, a correct explanation 
of pr puzzle is expected to describe, quantitatively, or 
semi-quantitatively, the properties of all measured decay 
modes. 

At last, a few words about Fock state component. This 
concept has been discussed in many works, such as by 
BK [41], CB [36], and FK [79]. Even from the theoreti- 
cal explanation of experiments other than et e™~ collision, 
the non-perturbative Fock component is also indicated to 
exist. For example, the analysis of charmonia photopro- 
duction amplitude implies [80] that the |c€) component 
is narrowly distributed in the transversal direction, while 
the cé pairs at larger separations may be part of higher 
Fock states which contain gluons and light quarks. 


II. EVALUATION OF RATIO OF y’ TO J/y 
DECAY 


In this section, we depict three approaches to estimate 
the ratio of y’ to J/w decay, one from theoretical analy- 
sis [1], the other two from experimental evaluation [81]. 


A. Theoretical Method 


Conventionally, the measured ratios of y’ to J/ẹ 
branching fractions for specific exclusive hadronic decays 
are compared with the naive prediction of pQCD, the so- 
called “12% rule”. In the framework of pQCD [1], the 
w particles are considered to be nonrelativistic bound 
states of a charm quark and its antiquark. Their decays 
into light hadrons are believed to be dominated by the 
annihilation of the cē pair into three gluons. In order 
to annihilate, the c and € must have a separation of or- 
der 1/m,, which is much smaller than the size of the 
charmonium state. Thus the annihilation amplitude for 
an S-wave state like J/w and w~’ must be proportional 
to the wave function at the origin, y(r = 0) [1]. The 


width for the decay into any specific final state h consist- 
ing of light hadrons is therefore proportional to |7(0)|?. 
The width for the decay into ete~ is also proportional 
to |w(0)|?. This leads to the simple prediction that the 
ratio of the branching fractions between y’ and J/w is 
given by Eq. (1). 

However, this naive prediction suffers from several ap- 
parent approximations. Higher order corrections, which 
may not even be small, are not included in this calcula- 
tion. For example, a first order correction to the branch- 
ing fraction of J/y — ete” could be 50% of the lowest 
term if as(my/y) ~ 0.2 [82] is taken. The relativistic 
effect is also ignored. Since the mass difference between 
J/w and y' is around 20% and < v?/c? >~ 0.24 for 
J/w, this correction may be at the same level as the low- 
est order contribution [82]. The inclusion of the finite 
size of the decay vertex will significantly reduce the ggg 
decay width of the J/w [83]. Moreover, the effect of non- 
perturbative dynamics is neglected, the size of which is 
hard to estimate. Therefore, people may question the 
validity of the 12 % rule as a serious benchmark for com- 
paring experimental data. 


B. Experiment Estimation (I) 


The first experimental estimation is based on the as- 
sumption that the decays of J/w and y’ in the lowest 
order of QCD are classified into hadronic decays (ggg), 
electromagnetic decays (y*), radiative decays into light 
hadrons (ygg), and transition to lower mass charmo- 
nium states (ccX) [72, 81]. Thus, using the relation 
B(ggg) + B(ygag) + B(y*) + B(X) = 1, one can derive 
B(ggg) + B(ygg) by subtracting B(y*) and B(X) from 
unity. 

The calculated values of B(y*) and B(céX), together 
with the values used to calculate them are summarized in 
Table III. As regards to w’, two final states yn(25) and 
h-(1+P,)+X with faint branching fractions are neglected 
in our calculation. By deducting the contributions B(7*) 
and B(X), we find that B(J/w — ggg) + B(J/w > 
799) = (73.3 £0.5)% and B(Y' — ggg) + BW > 799) = 
(18.9 + 1.3)%, then the ratio of them is 


_ BY" > 999 + 199) 
B(J/Y > 999 + 799) 


The above estimation is consistent with the previous 
ones [69, 81]. The relation between the decay rates of 
ggg and ygqg is readily calculated in pQCD to the first 
order as [82] 


Qg = (25.7 +1.7)%. (15) 


T(J/p— 799) _ 16 a as 
T(J/Y > ggg) 5 as(me) (1 ). 


Using as(m-) = 0.28, one can estimate the ratio to be 
0.062. A similar relation can be deduced for the w’ de- 
cays. So we obtain B(J/w — ggg) ~ (69.0 + 0.5)% and 
B(w! > ggg) ~ (17.8 + 1.2)%, while the “25.7% ratio” 


TABLE III: Experimental data on the branching fractions for 
J/ and w’ decays through virtual photon and to lower mass 
charmonium states used in this analysis. Most of the data are 
taken from PDG [27], except for B(J/w, Y’ — 7* — hadrons), 
which are calculated by the product R- B(J/y, Y > wtp), 
with R = 2.28+0.04 [84]. In estimating the errors of the sums, 
the correlations between the channels are considered [85]. 


Channel B(J/wW) 
y* — hadrons (13.50+0.30)% 


BW’) 
(1.66£0.18)% 


ete” (5.94+0.06)% (7.35+0.18) x 107% 
UTU (5.9340.06)% (7.3 +0.8)x 107? 
TtT — (2.8 +0.7)x 107° 
y >X (25.37£0.35)% — (3.41£0.27)% 
Me (1.340.4)% — (2.6£0.4)x 10-7 
ata Jj (31.8 +0.6)% 
non F/ op (16.46+0.35)% 
nJ/w (3.09 +0.08)% 
n? JJ (1.26+0.13)x 107° 
YXco (9.2+0.4)% 
YXe1 (8.740.4)% 
VXc2 (8.140.4)% 
COX (1.340.4)% (77.741.2)% 


in Eq. (15) stands well for both ggg and ygg. Although 
Qg is considerably enhanced relative to Qa in Eq. (1), it 
is fairly compatible with the ratios for the K+K- and 
K2 K? decay modes between y’ and J/w, which are 


Qrir- = (26.644.5)% (CLEO), 
Qrogo = (32.2 + 5.2)% (CLEO) , 


according to the recent results from CLEO and BES [5- 
7]. The relation in Eq. (15) was discussed in the literature 
as the hadronic excess in w’ decays [69, 81]. It implicates 
that while some modes are suppressed in 7)’ decays, the 
dominant part of y’ through ggg decays is enhanced rela- 
tive to the 12% rule prediction in the light of J/w decays. 


C. Experiment Estimation (II) 


The second approach for estimating Qa is to use 
the data on branching fractions for hadronic decays 
in final states containing pions, kaons, and protons 
that have already been measured for both the J/y 
and the Yy’. They are ntn, KEK, pp, mtx n°, 
pon, (ntr), 3(ata-), (rtr )n?, 3(rta—)n°, 
Antr), (KTK), K*K-ata-, K Knnn’, 
KEKA ntr), rtn pp, ppr*r n? and so forth. Us- 
ing the data compiled in Table IV, we have 


XC B(J/b > fi) = (14.46 + 0.80)% 


i=l 


and 
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SOBW > fi) = (1.70 £ 0.23)% . 


i=l 


It follows that 


Qs = Y BW > fi) 


TABLE IV: Branching fractions for 7)’ and J/ decays, and 
Qn values are also calculated. For y’ > ata, K+K- de- 
cays, the results are the weighted average of measurements 
from CLEO [5] and BES [86], while for J/7) —> nrn’, the 
result is the weighted average of measurements from BES [3] 
and BABAR [4]. Except for values with f, all others from 
PDG [27]. 


Final state Byjy(10~*) Byes (10) Qa (%) 
nTn 0.147 + 0.023 0.08 + 0.05f 5.63.5 
KK- 0.237 + 0.031 0.63 +0.07f 26.5 + 4.5 
pp 2.17 +0.08 2.65+0.22 12.2 1.1 
ata n? 21.2 + 1.01łf 1.68 +0.26 0.79 0.13 
ppr? 1.09+ 0.09 1.3340.17 12.2+ 1.9 
(a'r) 3.55 + 0.23 2.4+0.6 6.8+£1.8 
3(a' r) 4.3+ 0.4 3.5 + 2.0 8.1+4.7 
3(at a7)? 29+6 35+ 16 12.1+6.1 
(ntr )n° 33.7 + 2.6 26.6 + 2.9 7.91.1 
(rtnn?) 16.2 + 2.1 4514 27.8 + 9.4 
2(Kt K7) 0.78 +0.14 0.60+0.14 7.72.3 
KK- ntr” 7.2 + 2.3 T2405 10.0433 
KHK ntr’? 12.0+3.0 12.4+1.0 10.32.7 
KtK-2(nta-) 4.70.7 18 + 38.3 + 20.0 
pprt 1 6.0 £ 0.5 6.0 + 0.4 10.0 1.1 
pprtnr 7° 2.30.9 7.30.7 31.7 + 12.8 


Some remarks are in order here. First, we know that 
most of the multihadron final states in fact include sums 
of several two-body intermediate states, so the Qs is 
not the exact ratio of y’ to J/w inclusive hadronic de- 
cay rates, but the ratio on average of the exclusive de- 
cay channels as measured to date. In another word, 
Qs represents a mixed effect which may deviate notice- 
ably from the component Q-values. For example, the 
decay Y — ntr K+K can proceed through inter- 
mediate state K*(892)°K5(1430)° + c.c., whose Q = 
(2.9 + 1.3)%, is greatly suppressed [19] comparing with 
Q = (10.0 + 3.3)%. Second, we notice that the results 
obtained by two estimations vary considerably; further- 
more, by virtue of Table IV, it is obvious that many 
Q-values deivate from Q, significantly while on average, 
Qs is similar to pQCD Qna. It would be an intriguing 
problem that with more and more data and higher accu- 
rate measurements if Q, and Qg could approximate with 
each other and be consistent with Qa, or if Q, and Qg 
still deviate from each other prominently. 


IV. COMMENT ON THE pr PUZZLE 
A. Theory Aspect 


For all explanations involving pr puzzle, there is a ba- 
sic assumption that the non-relativistic potential model 
delineates the physics of charmonium decays to a good 
approximation. However, this assumption indeed re- 
quires examination in detail. First, as we have noticed 
in subsection II A, several corrections, should be added 
for the decay ratio between yY’ and J/w. 

Second, we see the effect of non-perturbative dynam- 
ics is neglected, which is crucial for charmonium decays. 
Actually from various explanations of pr puzzle, or more 
generally from the phenomenological explanation involv- 
ing pQCD, we note that certain non-perturbative effect 
or nonlinear effect must be incorporated one way or the 
other in order to recount the experimental data. But this 
kind of effect could hardly been included in the present 
non-relativistic potential model. 

Third, if the S-D wave mixing scenario holds as the 
solution of pr puzzle, then the matrix element of D wave 
to light hadrons would be very large, which can hardly be 
accommodated in the potential model, or any other cur- 
rent theory. It indicates that the current understanding 
of charmonium decays may not be complete. 

We envisage that new development of the theory 
should take into account the following features involving 
charmonium decays 


e the mass effect for different charmonium states, 
such as J/w, ~’, Y” and so on; 


e the non-perturbative or non-linear effect on the res- 
onance decays; 


e reasonable description for the known features of 
charmonium spectroscopy; 


e quantitative consistence (with reasonable high ac- 
curacy) with the existing experimental measure- 
ments. 


Herein it is also important to distinguish the quarko- 
nium states (theoretical states) and mass eigenstates 
(physical states) [51, 53]. From a fairly theoretical point 
of view, if both states, say, J/y and w’ did dominantly 
annihilate into three gluons, they would mix and could 
thus not be the putative quarkonium states. In a non- 
relativistic potential model, for example, the w’ is simply 
a radial excitation of the J/y. This is a well defined pic- 
ture in which J/w and y’ are orthogonal states. If the 
annihilation into three gluons could be treated as a “per- 
turbation” to the non-relativistic potential, then clearly 
the unperturbed states would mix and rearrange them- 
selves into orthogonal mass eigenstates. The QCD dy- 
namics may be such that the physical states, presumably 
mixtures of the theoretical cé@ states (n?5+!L 7), are so 
built up that one of them strongly annihilates into three 
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perturbative gluons while the other does not. Mixing of 
the 13S, and 23S; states via three perturbative gluons 
has little effect on the charmonium mass spectrum, but 
may be crucial for the decay pattern. 


B. Phenomenology Aspect 


One may remember that at the early stage of the dis- 
covery of a narrow state J/w, the cé system was hailed as 
the Hydrogen atom of Quantum Chromodynamics, with 
the implied hope that the study of the newly discov- 
ered system could shed as much light on the dynamics of 
quark-antiquark interactions as the study of the Hydro- 
gen atom had on Quantum Electrodynamics. But one 
may also notice the historical fact, even before Bohr’s 
theory, Balmer series had been discovered for long, and 
the famous Rydberg formula had also been proposed, 
which laid a solid foundation for further theoretical im- 
provement. If we are conscious of the more complicat- 
edness of charmonia system comparing with Hydrogen 
atom, we may prepare for more hard and meticulous 
works. As a first step, it is necessary to develop a reli- 
able and extensively applicable phenomenological model 
(PM). 

The advantage of PM lies in that a well-defined 
PM contains few experimentally determined parameters 
which have clear physical meaning; moreover, with only 
few parameters determined from experiment, PM could 
produce concrete results which can be directly confirmed 
or falsified by experiment and may guid experimental 
searches. Such a model has a good relation with elemen- 
tary principle of the theory, and if correct, can be used 
for further theoretical refinement. This point is note- 
worthy for the time being, since the general QCD can 
hardly provide solutions for special problems; conversely, 
we have to establish certain effective empirical model to 
advance our understanding for generic QCD principle. 

Here we would like to mention few ideas of PM, which 
have or intend to provide a general description for char- 
monium decay. 


1. Mixing Model 


As we have noticed in subsection II C4, whatever CI 
model or FK model, they could yield definite predic- 
tions for experimental test, and therefore provide clues 
for further development. In addition, according to Haber 
and Perrier’s parameterization scheme [77], we could get 
the decay rate relations between different channels, and 
the proportions of different interaction amplitudes, all of 
which are useful information to understand the dynamics 
of charmonium decays. 

The special feature of Rosner’s S- and D-wave mixing 
scheme [56] is that it is simple, and it works both for 
suppressed and enhanced decay modes, and moreover, it 


connects the J/q, y’ and w” decays together, and give 
predictions with little uncertainty. 


2. Effective FSI 


Miller once discussed the connection between the 
strong-coupling approximation to quantum chromody- 
namics and nuclear properties observed at low and 
medium energies and momentum transfer [87]. He sug- 
gested that the strong-coupling (corresponds to long- 
range interaction) version of QCD does reproduce the 
salient feature of the meson-baryon picture of low mo- 
mentum transfer in nuclear physics. The derivation of 
nuclear physics indicates that quark aspects of ordinary 
nuclei are hidden in the hadronic degrees of freedom. Re- 
versely, one may imagine that a quasi-meson or quasi- 
baryon structure could be formed in non-perturbative 
hadronization process. The meson-baryon picture in nu- 
clear physics could be utilized, with some modification, as 
an effective FSI theory, or molecular-model theory. The 
recent molecular interpretations [88, 89] for the newly 
found state Y (4260) could be treated as such an effective 
FSI theory. 

Furthermore, the quasi-particles could even be real 
particles, and then with the residual strong force between 
the quarks inside the quasi-particles, some multiplets can 
be formed as the quasi-particle and anti-quasi-particle 
bound states. This idea was first put forth to explain a 
lots of low mass baryon-antibayron enhancements newly 
found by CLEO, Belle and BaBar collaborations [90]. 


3.  Glueball and Hybrid 


We have seen that the J/w-glueball mixing scheme is 
the first explanation proposed to explain the pa puzzle. 
With the implication from lattice calculation [91], Suzuki 
once proposed a glueball-w’ mixing scheme to explain the 
excess hadronic decays at the y’ [92]. In fact, gluodynam- 
ics is always a tantalizing domain for theorists. Recently, 
some lattice evaluations suggest [93, 94] that the masses 
of some hybrid states could be low enough to be degen- 
erated with charmonium states, such as J/q and/or Xes. 
Unlike Glueball, these ceg hybrid states could be very 
narrow, so it may be very stimulating and noteworthy to 
search for and confirm these kinds of states. 


C. Experiment Aspect 


Physics is a science of experiment. Physical facts are 
bases for theoretical development and also criteria for 
checking theoretical hypotheses. It is a prominent fact 
that the early analyses on pa puzzle based on meager ex- 
perimental data often lead to unsatisfactory, sometimes 
premature, inferences, which were washed out easily by 
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later accurate data. In fact, the current knowledge con- 
cerning the Y’ decays from experiments is still rather lim- 
ited, even summing all charmless channels presented in 
PDG [27], the total branching fraction is less 2%. Such 
situation prevents us from laying down a solid founda- 
tion for elementary dynamics exploration. However, the 
estimation of Q value, as discussed in Section III, affords 
us some clues concerning the exploration of charmonium 
decay dynamics. Since many suppressed channels have 
been found, especially those such as pr which is greatly 
suppressed in w’ decay, and if the Qpa really represents 
the averaged value of inclusive hadronic decay, the esti- 
mation of Q, indicates that either lots of enhanced de- 
cays are not discovered, or some particular decays only 
present in Y’, or both cases exist. Therefore, systematic 
experimental study of w’ decays is anxiously awaited. 

Moreover, study should be carried out not only on the 
ratios of Y’ to J/w decays, but also on the other ratios 
such as those between n, and ne [79], and those between 
wp” and J/[50], and/or many other ratios between differ- 
ent resonances for the same channel or between different 
channels about the same resonance [54]. All these studies 
will shed light on the understanding of charmonium de- 
cays. From our point of view, the progress would be more 
likely obtained from the analysis of the new experimen- 
tal results, rather than from the inspiration of general 
theoretical principle. 


V. SUMMARY 


In this paper, we present a general review on the study 
of the pr puzzle. Although there is still no satisfactory 
explanation for all existing experiment results, some en- 
lightenment ideas have been put forth. In addition, we 
also discuss three methods of estimating the ratio of the 
branching fractions in J/~ and y’ decays. In the light of 
the present theoretical and experimental status, we argue 
that it is important to explore the potential models from 
a new point of view, it is necessary to search for and/or 
construct an effective phenomenological model, and it is 
especially crucial to perform systematic measurements of 
various charmonium decays. 

As we know physics in the charm energy region is in 
the boundary domain between perturbative and nonper- 
turbative QCD. Recently the observed hadronic decays of 
charmonium may give new challenges to the present theo- 
retical understanding of the decay mechanisms. Many of 
the new charmonium states observed by Belle and BaBar 
which can hardly be accommodated by potential models 
have led to new theoretical speculations about charmo- 
nium spectroscopy and novel production mechanism [95]. 

Experimentally, the expected large data sample from 
CLEOc in the near future, and even larger sample from 
forthcoming detector BESIII will open for us a new era 
for charmonium dynamics study, even though we may not 
obtain a throughly revolutionary theory, we could acquire 
more profound understanding for the existing theory, at 


the same time we could expect a brand-new comprehen- 
sion for the pr puzzle. 
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